Insulating gap in FeO: correlations and covalency 
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We report calculations of the electronic structure of FeO in the LDA and LDA+U approximation 
with and without rhombohedral distortion. In both cases LDA renders an antiferromagnetic metal, 
and LDA+U opens a Hubbard gap. However, the character of the gap is qualitatively different in 
the two structure, and the difference can be traced down to underlying LDA bnad structure. An 
analysis of the calculations gives a new insight on the origin of the insulating gap in 3d monoxides 
and on the role of the k-dependency of U, missing in the contemporary LDA+U method. 
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It is a well-knovi'n fact that conventional band struc- 
ture calculations incorrectly give metallic ground state 
for the intermediate 3d transition metal monoxides, CoO 
and FeO. Before the high- Tc cuprates entered the scene, 
this had been often considered as the most notable failure 
of the Local Density Approximation (LDA). In the last 
decade a number of extensions of the density-functional 
theory were suggested, which, in different manners, led to 
insulating ground states for this compounds. Most suc- 
cessful were varicws flavors ef the Self-Interaction Cor- 
rected (SIC) LDAEI, LDA+UEI, and the Orbital Polariza- 
tion Correction in a crystal field basisH. Interestingly, 
apart from the gap itself, LDA appears to do as good, 
and sometimes a better job, than these sophisticated 
extensions, especially when Generalized Gradient Cor- 
rection to the conventional LDA is taken into account. 
Structural properties are reproduced- very well, includ- 
ing rhombohedral distortion in FeOcI, and its increase 
with pressureO. Moreover, angular resolved photoemis- 
sion renders the bands more similar to the LDA bands, 
than to those in other calculatipnal scheme (except for 
the narrow range near the gap)Q. The magnitude of the 
magnetic moment, which is often believed to be the first 
indication of an LDA failure, is nearly exact for the spin 
moment in CoO (2.4 (Ib) and in FeO (3.5 /x^, assuming 
no orbital moment). On the contrary, the gap- improving 
calculation tend to overestimate the Fe moment, espe- 
cially wherL-the orbital moment is included (i.e., there is 
a tendencytl to underestimated the crystal field quench- 
ing of the orbital moment). It seems, then, that the main 
problem in the LDA is purely spectroscopical (quite in 
the spirit of the density functional theory), namely non- 
existence of the gap. It is worth noting that an important 
common feature of the LDA and the corrected schemes 
mentioned above is substantial width of the metal ri- 
bands (for instance, in LDA the width of Fe t2g band 
at normal pressure is about 1.4 eV). This should be con- 
trasted with the popular analysis of the electronic struc- 
ture of 3(i-oxides in terms of separate levels of a width 
less than Jan- Teller energy and spin-orbital couplingE. 



The fact that LDA gaps are always too small is well 
understood. Mathematically it appears as the density- 
derivative discontinuity of the Kohn-Sham potential in 
the DFT. The physics of this discontinuity may be dif- 
ferent, but for transition metal oxides (NiO, MnO) ib 
is usually associated with the Mott-Hubbard repulsionEl. 
In view of this, it has always been much more disturb- 
ing to have wrongly a metallic behavior in LDA calcula- 
tion, thap just to have a wrong gap. As formulated by 
Normand, "one would like to obtain a gap at the level 
of a density functional calculation (no matter how small) 
so as to define the Mott-Hubbard correction in an unam- 
biguous fashion" . Moreover, even in a case when LDA 
does not give a gap, but gives reasonable band structure 
except for the immediate vicinity of the Fermi level, and 
correctly describes delicate features of the ground state, 
like magnetoelastic interactions, it is desirable to have a 
correction scheme which does not destroy the LDA bands 
completely, but rather corrects them in a systematic 
manner. Unfortunately none of the schemes above acts in 
such a way. In this paper we shalli-analyze the results of 
the rotationally invariant LDA+UlI calculations for FeO 
in more details than is usually done, and compare them 
with the standard LDA calculations, paying particular 
attention to the process of the LDA+U gap opening in 
cubic and rhombohedral structure. We will see that the 
ground state in the LDA+U approach is intimately re- 
lated with the underlying LDA band structure, although 
LDA+U cannot fully account for the bands hybridization 
effects, which seem to be quite important here. We will 
argue that none of the existing "corrected LDA" schemes 
(nor the straight LDA) correctly describes the insulat- 
ing ground state in FeO and similar compounds. On 
the other hand, non-local schemes similar to the GW 
approximation may be able to provide a qualitatively 
correct description. It is worth noting that none of the 
existing "extended-LDA" calculations has taken into ac- 
count such an important factor as the distortion from 
the ideal cubic NaCl-type structure, which is associated 
in 3(i-monoxides with the onset of magnetic ordering. 



1 



Electronic structure of an isolated Fe"'""'" ion in a cubic 
field is [described in the high-spin state by the following 
schemeEl: The spin-up d-states are all filled, and sepa- 
rated from the spin-down states by the exchange splitting 
Eex- The partially occupied spin-down states are split 
by the crystal field, so that e{eg) — e{t2g) = A ^ E^x- 
There is one electron in the t2gi state, which is triply de- 
generate. Antiferromagnetic FeO has close-packed (111) 
planes of Fe ions with the same spin, which alternate with 
the similar planes with the opposite spin, thus lowering 
the symmetry to rhombohedral one. In the rhombohe- 
dral field, there is one linear combination of the t2g states, 
namely Aig — {xy + yz + zx)/^/3, which can also be writ- 
ten as — r^, where Z is parallel to [111] (here and be- 
low we shall use lower case symbols for the states classifi- 
cation in cubic symmetry and upper case symbols for the 
rhombohedral nomenclature). The four other states have 
the same symmetry Eg (in the rhombohedral nomencla- 
ture), but when the deviation from the cubic symmetry 
is small one can speak about the upper two levels. Eg , 

and the lower two levels, Eg, which are close to the Aig 
also originate from the cubic t2g states, and are separated 
from Eg by approximately A. In the LDA calculation, 
when the crystal symmetry is still cubic, although the 
magnetic ordering is rhombohedral, the splitting between 
Eg and Aig is small, much smaller than their bandwidths, 
so that they merge into one band (corresponding to the 
cubic t2g band), which is necessarily metallic. Note that 
although the magnetism in this system appears due to 
the indirect exchange and is determined by the 0-Fe pd 
interaction, the width of the t2g band is mainly due to 
the Fe-Fe dda overlap. Fe-0 hopping for this band is 
mainly pdn and weak. The situation in CoO with two 
t2g electrons is very similar. 

There were several sujCCjessful attempts to obtain an 
insulating state in FeOETO. Interestingly, all these ap- 
proaches give insulating gap in a fair agreement with 
the experirrient, but all because of different reasons. SIC 
calculationsld favor d— bands in pure orbital states, that 
is, undermining the role of crystal field and quenching 
of orbital moment. The effect on the occupied d-states 
is thus extremely strong (of the order of 1 Ry) and all 
four oxides come out as pure charge transfer insulators. 
This is in contradiction with the general experimental in- 
dication that the character of the band gap changes from 
predominantly Mott-Hubbard to predominantly charge- 
transfer when going from MnO to NiO. Large orbital 
moment is obtained for FeO, in contradiction with the 
experiment (an argument is usually made that the ex- 
perimental number may be incorrect because poor sam- 
ple quality). It is also worth noting that the SIC for- 
malism was initially invented as a remedy of the LDA 
in the direction of the exact, self-interaction free density 
functional. It is hard to imagine, however, that the ex- 
act density functional theory with its orbital-independent 
one-electron potential can unquench the orbital moiHcnt. 

The orbital-corrected functional used by NormanB has 



a similar problem. The correct many-body solution fojL 
an isolated ion with unfilled t2g shell in a cubic fieldB 
has an energy contribution proportional to the total an- 
gular momentum. This term, however, does not favor 
a specific direction of the momentum, that is, does not 
include the projections. The only interaction which does 
unquench the orbital momentum is spin-orbit, and it is 
relatively weak. Norman's correction has the same func.- 
tional form as the exact quantum chemical expressionQ, 
but substitutes the total momentum by its projection on 
the quantization axis. This orbital moment projection 
dependent term acts in a way similar to the spin-orbit 
coupling, but with a much larger magnitude (of the or- 
der of exchange splitting J). Thus in FeO the occupied 
spin-down d-band is too close to a pure m = 1 state, that 
is, xz + iyz. 
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FIG. 1. (a) LDA band structure without a rhombohedral 
distortion. The bars show the relative Aig character of the 
correspoding states; (b) LDA+U bands for the same crystal 
structure. Note the bottom of the Fe(s) band at the F point. 

The corresponding band structure is shown on Fig||, 
where the LDA+U, equally successful in predicting a 
band gap, forms a completely different insulating state. 
There is no mechanism for unquenching the orbital mo- 
menL, apart from the spin-orbit coupling, neglected in 
Ref.H and other LDA-fU calculations. LDA-I-U starts 
from the straight LDA bands, where the tig bands are 
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strongly mixed, the Aig (in rhombohedral notations) or- 
bital being slightly more occupied than each of the two 
Eg orbitals. When formulated in rotationally invariant 
formcl, LDA-t-U tends to apply the positive U correction 
to less occupied Eg orbitals, making them less and less 
filled (potential U correction for more filled Aig-orbital is 
negative) , and eventually splits off the Aig band, forming 
a gap between it and the Eg bands. Sufhciently large U 
(we used the empirical value of C/=5.1 eVliJ, which giv£a 
a good value for the gap; constrain LDA calculationstiJ 
yield a somewhat larger U—6.8 eV) pushes the occupied 
Aig band down close to, and for U—6.8 eV right into, 
the 0(p)-bands manifold, and the unoccupied Eg bands 
up above the bottom of the Fe(s) band. This seems to 
be in qualitative agreement with the photoemission ex- 
periments, which show the top of-the valence band to be 
of the mixed 02p-Fe3d characterE£l, and with the optical 
experiments, which indicate a weak absorption between 
0.5 and 2.0 eV, assigned to the (pd) — s transitions, and 
strong absorption edge at 2.4 eV due to transitions into 
the Fe{d) bandB. In our calculations, the minimal gap 
opens between the 0{p)—Fc{d) band and the Fe(s) band 
at 1.3 eV and transitions into the Fe{d) band start at 2.2 
eV. 

LDA and the LDA-I-U bands refer to the cubic-jFeO 
(here and below we show results of the ASA LMTOEl cal- 
culations with the unit cell volume 256 bohr'^/FeO, close 
to the LDA equilibrium volumoj). Since only 1 orbital, 
Aig, is occupied (see Table ||), there is no orbital mo- 
ment. One can expect that including weak (a few mRy) 
spin-orbit coupling will create a small orbital moment, 
but hardly one comparable with the spin magnetization. 
It is interesting that the way the gap opens in LDA+U in 
the cubic structure is very "LDA-like" : The only way to 
open a gap in an effective one-electron approximation is 
to split the t2g band by (magnetic, in this case) rhombo- 
hedral symmetry and to occupy the Aig orbital for FeO 
and the two Eg orbital in CoO. In fact, this is exactly 
what happened in the LDA calculations of Dufek et allJ, 
who used an exotic LDA functional, which gave poor to- 
tal energies but did open gaps in both compounds. 

Overall this seems to be a physically satisfactory de- 
scription of the gap opening. However, this is not as 
straightforward as one may think. To show the prob- 
lem, let us compare the calculations in the cubic struc- 
ture with those in the rhombohedrally distorted struc- 
ture (Fig. ^ we used the LDA equilibrium distortion, 
calculated in Refs. |^; experimental distortion is smaller, 
but increases with the pressure). To understand the re- 
sult, one should keep in mind that dispersion of the t2g 
bands is mainly due to direct dda hopping between the 
like spins, which can be easily verified in the LMTO-TB 
method by removing oxygen orbitals from the basis set, 
or by looking directly on the corresponding elements of 
the LMTO-TB Hamiltonian. We observe that the dis- 
tortion increases the t2g bandwidth (due to decreased 
distance between the like-spin ions), but mainly at the 



expense of the E^-like bands. The dispersion of the Aig 
band decreases instead, thus leading to a decreased occu- 
pancy of the Aig-like state (because the whole t2g mani- 
fold is less than half filled) . One can see this by compar- 
ing the occupancy matrices in the cubic (Table ^ and in 
the distorted (Table III) structure on the first iteration, 
i. e., before the effect of U. While in the cubic structure 
the Aig state was the mostly occupied one (with the two 
Eg states close next), in the distorted structure the Eg 
states are twice more occupied than the Aig one. 
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FIG. 2. The same as Fig.l, but with the rhombohedral 
distortion of 6%. 

The LDA-I-U code we used is rotationally invarianti, 
which means that the program can apply the U correc- 
tions to the states which are arbitrary linear combina- 
tions of the d-orbitals with various to. The coefficients 
of this combination are chosen self-consistently so that 
they optimize the total energy including Coulomb one. 
In the case of the rhombohedrally distorted FeO it is en- 
ergetically unfavorable to occupy the Aig orbital, but not 
because of the Hubbard energy, but because of the one- 
electron energy: the same reason why LDA deoccupied 
this orbital after distortion. How does LDA+U handles 
this situation during the self-consistency cycle? It is left 
with the only possibility to make a linear combination 
of the Eg orbitals with a possible admixture of the Eg 
states, as much as the crystal field allows that, to make 



3 



it a separate band and to open a gap between this band 
and the rest of the spin-minority c?-bands (Table [V). 

One may feel it unphysical that the method converges 
to qualitatively different ground states with and without 
distortion. This fact suggests that probably the correct 
solution should be intermediate between the two LDA+U 
stationary points, a completely full and a completely 
empty Aig state. Indeed, a closer look at the LDA band 
structure of the Fig.||a shows that along the whole F-Z 
(111) direction the Aig band lies entirely below the Fermi 
level. In the other parts of the Brillouin zone (points L 
and F) it is completely above the Fermi level. Thus, a 
U matrix acting on an unoccupied Aig orbital has to 
overcome the natural LDA tendency along the F-Z line, 
although it is in accord with the LDA tendencies in the 
rest of the zone. 

From this consideration one can see that the problem 
with the LDA+U is that the [/-matrix is k-independent: 



it is a matrix in the orbital space, U„ 



-2 < TO < 2, 



and this matrix does not know whether the particular d- 
electron is in a bonding or antibonding interaction with 
the neigboring sites. While this may be a good approx- 
imation for extremely localized electrons, in such mate- 
rials as FeO one cannot neglect the fact that d-orbitals 
extend into neighboring sites. In the LMTO scheme the 
tails of each orbital, penetrating into an atomic sphere of 
a neigboring site, is re-expanded and, depending on the 
LMTO flavor ("representation"), may or may not ap.- 
pear in the density matrix at the same site as its headll3. 
In the conventional LDA-I-U scheme, as well as in all 
other correction schemes discussed above, only the heads 
of the LMT orbitals are subjects to a correction. On 
the other hand, a more realistic correction would also ap- 
ply to that part of a d-orbital that penetrates into the 
neighboring spheres. In other words, a physically correct 
scheme for implemented such a correction should be es- 
sentially non-local. Hartsee-Fock like schemes, similar to 
the GW approximationim, should be able to reproduce 
the correct physics. For instance, one can expect that 
the occupied band that appears to be purely Aig in the 
cubic case and purely Eg in the distorted case, would 
have a mixed character, being more ^i^-like close to the 
F-Z line. Correspondingly, the effective [/-correction ma- 
trix would be different when applied at different k-points. 
Another interesting alternative is offered by non-Kohn- 
Sham versions of the Density Functional Theory, where 
a non-local (screened Hartree-Fock like) contribution to 
the total energy is singled out together with the non- 
interacting kinetic energy, and the rest is treated in an 
LDA. Such schemes retain the good accuracy of the total 
energy calculations in the LDA, and improve the excita- 
tion energies as well. It is worth noting that underestima- 
tion of the hybridization tendencies inside the tig band in 
LDA-I-U and other schemes is likely to lead to a wrong en- 
ergetics and spoil the agreement of structural properties, 
calculated in LDAD, with the experiment. More detailed 
discussion of structural and magi^toelastic properties of 
FeO will be published elsewhereta. 



To summarize, we report conventional LDA calcula- 
tions and rotationally-invariant LDA-I-U calculations for 
antiferromagnetic FeO, both in the cubic and in a rhom- 
bohedrally distorted structure. In both cases LDA+U 
opens gaps, but these gaps are of completely different 
character in the two structures, and also different from 
the gaps appearing in other "corrected-LDA" schemes, 
which, in turn, differ by their physical nature from each 
other. We believe that this is a consequence of the lo- 
cal character of the Mott-Hubbard correction in the con- 
ventional LDA+U method, and neither straight LDA, 
nor "corrected-LDA" methods provide a proper physi- 
cal description of the gap formation. We argue that the 
ultimate method must take into account the non-local 
character of the Coulomb repulsion. 
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TABLE I. Orbitals and their occupations without rhombo- 
hedral distortion, for U = 0. Coordinates system corresponds 
to the rhombohedral symmetry: z-axis is perpendicular to 
the ferromagnetic Fe planes, y-axis points towards the nearest 
neighbor Fe with the like spin. The character of the orbitals 
is, in the order of the table, 2 x Eg, 2 x Eg, Aig 



occupation 


xy 


yz 


3z^ - 1 


xz 




.17 


-.03 


.03 


.00 


-.65 


.76 


.17 


.75 


-.66 


.00 


-.03 


.03 


.31 


.60 


.69 


.00 


.30 


.26 


.33 


-.26 


-.30 


.00 


.69 


.60 


.39 


.00 


.00 


1.00 


.00 


.00 



TABLE II. Orbitals and their occupations without rhom- 
bohedral distortion, for U = 5.1. The character of the orbitals 

is 2 X Eg, 2 X Eg, Aig 



occupation 


xy 


yz 


3z^ - 1 


xz 


2 2 

X -y 


.04 


-.10 


-.10 


.00 


.68 


.72 


.04 


.72 


.68 


.00 


.10 


.10 


.10 


.53 


-.57 


.00 


-.46 


.43 


.10 


-.43 


.46 


.00 


-.57 


.53 


.90 


.00 


.00 


1.00 


.00 


.00 



TABLE III. Orbitals and their occupations with 6% 
rhombohedral distortion, for U = 0. The characters are 
2 X E'g, Aig, 2 X E'J 



occupation 






Wz- - 1 






0.17 


-.55 


.84 


.00 


-.02 


.01 


0.17 


-.01 


.02 


.00 


.83 


-.56 


0.19 


.00 


.00 


1.0 


.00 


.00 


0.38 


.84 


.55 


.00 


.00 


-.00 


0.39 


-.01 


-.01 


.01 


.56 


.83 



TABLE IV. Orbitals and their occupations with 6% rhom- 
bohedral distortion, for U = 5.1. The characters are 

^ISl -^S'j 2 X Eg, Aig, Eg 



occupation 


xy 


yz 


3«^ - 1 


xz 


2 2 

X -y 


.04 


-.02 


.00 


.99 


.14 


.04 


.04 


.96 


.28 


.02 


.01 


.01 


.10 


-.28 


.96 


.00 


.02 


-.01 


.10 


.01 


-.02 


-.09 


.83 


-.56 


.86 


-.01 


.00 


-.12 


.55 


.83 
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